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Abstract 

The long time behavior of an evaporating black hole presents a chal- 
lenge to theoretical physics and touches relevant conceptual issues of 
quantum gravity, such as the information paradox. There are basically 
two strategies: top-down, i.e., to construct first a full quantum theory of 
gravity and to discuss black hole evaporation as a particular application 
thereof, and bottom- up, i.e., to sidestep the difficulties inherent to the 
former approach by invoking "reasonable" ad-hoc assumptions. 

Exploiting the fact that the Schwarzschild black hole can be described 
by means of an effective theory in 2D, a particular dilaton gravity model, 
the latter route is pursued. A crucial technical ingredient is Izawa's result 
on consistent deformations of 2D BF theory, while the most relevant phys- 
ical assumption is boundedness of the asymptotic matter flux during the 
whole evaporation process. Together with technical assumptions which 
can be relaxed, the dynamics of the evaporating black hole is described 
by means of consistent deformations of the underlying gauge symmetries 
with only one important deformation parameter. An attractor solution, 
the endpoint of the evaporation process, is found. Its metric is flat. How- 
ever, the behavior of the dilaton field (which corresponds to the surface 
area) is nontrivial: it is argued that during the final flicker a first order 
phase transition occurs from a linear to a constant dilaton vacuum. Con- 
sequently, a shock wave is emitted as a final "thunderbolt" with a total 
energy of a fraction of the Planck mass. Relations to ultrarelativistic 
boosts are pointed out. Another fraction of the Planck mass may reside 
in a cold remnant. 

The physical discussion addresses the life time, the specific heat, the 
Carter-Penrose diagram, the information paradox and cosmological impli- 
cations. The phenomenon of "dilaton evaporation" to a constant dilaton 
vacuum might be of relevance also for higherdimensional scalar tensor 
theories. 
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1 INTRODUCTION 



1 Introduction 

For more than three decades the physics of black holes (BHs) has attracted 
the attention of an increasing number of relativists, astrophysicists, string the- 
oreticians and elementary particle physicists (for a recent textbook cf. e.g. [1]). 
One of the most spectacular theoretical predictions is the Hawking effect: due 
to quantum creation of particles a BH may evaporate with ensuing radiation 
having a thermal distribution [2] . 

However, there are well-known problems with the long-time evolution of 
BHs (e.g. in the form of the information paradox — ef. e.g. [3])). Although semi- 
classical models have established interesting insights, especially in the context of 
twodimensional (henceforth 2D) dilaton gravity (cf. e.g. [4] with particular focus 
on the model by Russo, Susskind and Thorlacius) the main obstacle in such a 
discussion stems from the fact that a comprehensive theory of quantum gravity 
does not yet exist. Moreover, already simple classical systems of gravity with 
matter are not exactly soluble which, in general, prohibits a nonperturbative 
treatment of quantum backreactions. 

The present approach tries to circumvent these difficulties at the cost of 
several ad-hoc assumptions, some of which could be called "natural" . The main 
idea is that one allows for a continuous family of models labelled by one or 
more deformation parameters. The word "deformation" will be given a precise 
meaning when the assumptions are specified, but roughly speaking one has a 
family of actions labelled by these parameters. The difhcult task is to establish 
an evolution equation for the latter. Dynamics then will dictate the evolution of 
these parameters. If an attractor solution exists dynamical deformations imply 
the transition of the original model, e.g. the Schwarzschild black hole (SS BH), 
to a specific remnant geometry, namely the attractor. Thus, rather than trying 
to find one model describing an evaporating BH from the early stages until the 
very end we have a family of models at our disposal and at each instant one 
family member provides the most adequate description of geometry. Provided 
the notion of light-like infinity {./) makes sense — this will be one of the ad hoc 
assumptions — one can imagine an asymptotic observer who patches all these 
models together along u = const, lines, where u is the retarded time. 

To be more concrete, the assumptions chosen in the present paper are 1. 
the restriction to consistent deformations in the technical sense of Barnich and 
Henneaux (roughly speaking, neither gauge degrees of freedom nor physical 
degrees of freedom may appear or disappear, but the gauge symmetries may 
be deformed, e.g. from an abelian gauge symmetry to a nonabelian one), 2. the 
assumption of asymptotic flatness of each spacetime solving the equations of 
motion of the deformed action during the whole process of evaporation, 3. the 
absence of nonextremal horizons to avoid extremal endpoints of the evolution 
(actually, a further technical assumption will be imposed which, however, can 
be relaxed), 4. an observer at spatial infinity (z*^) is assumed who measures a 
constant Hawking flux, i.e. the increase of energy loss of the black hole due to the 
Hawking effect is compensated by a corresponding boost of the observer. It is 
emphasized that this is merely a technical trick to allow for an easier derivation 
of the attractor solution. Of course, consequently practically nothing can be 
said about what happens at a given instant during the evaporation in terms 
of "physical time" — for this information not only the existence of such a boost 
but the precise value of the boost parameter would be needed. But since our 
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task is much more humble — we are only interested in the specific form of the 
remnant geometry and not a thorough description of the complete dynamics — 
these assumptions will be sufficient to provide a unique answer for the remnant 
geometry. 

The paper is organized as follows: the four assumptions relevant for the SS 
BH are presented in sect. |21 As a result due to the evaporation process the 
SS BH is deformed towards an effectively 2 + 1 dimensional model — in more 
picturesque (albeit slightly misleading) terms the celestial sphere evaporates 
into a celestial circle. A physical discussion of the results (e.g. the life time of 
a SS BH, its line element, the relevance of Planck scale contributions, dilaton 
evaporation, consequences for the information paradox and cosmological impli- 
cations) in sect. Eland possible relaxations of the assumptions (sect.QJ conclude 
this work. Complementary material can be found in the appendices: appendix 
^recapitulates deformations of generic dilaton gravity; a possible generalization 
which relaxes the technical assumptions is discussed in appendix A simple 
toy model from classical mechanics which mimics some of the main features of 
BH evaporation as described in the present work can be found in appendix ICl 

2 Evaporation of a Schwarzschild black hole 

The SS BH is the simplest of phenomenologically relevant BH models. Once the 
long time behavior of this case is understood it appears that generalizations to 
Reissner-Nordstrom or Kerr-Newman are comparably easy. More will be said 
on this in sect. 01 for now the focus will be on the SS BH. The formation process 
is completely neglected in these considerations. It is assumed that space-time 
has "settled down" after the collapse of matter and an approximately stationary 
situation is encountered. 

A "natural" working hypothesis, namely that no deformation occurs and 
spacetime remains SS for all times, leads to the usual prediction of an explosive 
evaporation during the final stage, being the result of applying the semi-classical 
approximation beyond its limit of validity, because eventually quantum backre- 
actions will deform geometry appreciably. 

It is noteworthy that no assumption is made regarding the type of matter 
coupled to gravity. It is just taken for granted that some matter degrees of 
freedom supporting the Hawking radiation actually do exist. Specifying matter 
explicitly has the disadvantage that one would have to find an exact solution of 
the coupled matter-gravity system in order to study the long time evolution of 
an evaporating BH. The present approach circumvents these difficulties at the 
price of several ad hoc assumptions.^ 

^It should be pointed out that the idea to sidestep the difficulties of quantum gravity and 
nevertheless obtain relevant predictions clearly is not new and has been pursued, for instance, 
in the framework of Doubly Special Relativity (for a recent review and more references cf. 
e.g. [7]) or in the context of /t-deformations [8]. Recent generalizations of some of these 
results to "gravity with an invariant energy scale" exist [9], but one has to be careful with 
the definition of the line-element [10,11]. A different philosophy which appears to be closer 
to the one of the present work consists of taking the Hilbert-Einstein action without further 
structure as starting point [12], The key ingredient in that reference is a running Newton 
constant. The causal structure is like that of a Reissner-Nordstrom BH because an inner 
horizon develops. The outer horizon decreases, like in the present approach; there seems to 
be a regular dS core (as suggested by [13]). However, as a careful analysis in their work shows 
nonetheless a singularity at r = is present, because the Newton constant scales differently 
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2.1 The four assumptions 

In many instances — take for example the singularity theorems in General Re- 
lativity — it helped to distinguish between the principal part of assumptions 
(e.g. one has to specify an energy condition) and their specific form (e.g. the 
strong energy condition). Similarly, one can relax or generalize the assumptions 
presented below. Some of these generalizations are addressed in sect. 0] and 
appendix^ On a first reading it might be useful to consult appendix lUl at this 
point. 

For each of the four assumptions first its realization in the toy example in 
appendix |0 will be mentioned, then the general form will be presented and 
finally its specification relevant to the core part of this work together with its 
implementation. 

2.1.1 First assumption: Consistent deformations 

Toy example: Due to interactions with the unknown "matter" part the geo- 
metric system (|C.l|l can be deformed by consistent deformations. 

General form: Due to quantum gravity effects during the long-time evolution 
of a BH, the gauge symmetries of gravity may be deformed in a specific way im- 
posed by hand (e.g. allow no deformations at all or only a certain class thereof). 
In principle, there can be an arbitrary amount of deformation parameters, but 
of course it is advantageous to keep its number as low as possible. 

Specific form: The BH evaporation due to the Hawking effect (plus even- 
tual non-perturbative effects from quantum backreactions) does not change the 
number of gauge degrees of freedom or physical degrees of freedom — i.e. the 
evaporation induces only consistent deformations^ of the underlying model (in 
the sense of Barnich and Henneaux [5]). Since spherical reduction and consis- 
tent deformations will not commute in general, it should be pointed out that 
consistent deformations are applied solely to the reduced theory. Thus, we are 
not dealing with the SS BH in the framework of Einstein gravity in = 4, but 
with its (classically equivalent) counterpart as 2D dilaton gravity. Obviously, 
deformations exciting nonspherical modes cannot be treated in such a simple 
manner. This situation resembles the so-called "spherical reduction anomaly" : 
spherical reduction and renormalization do not commute [18]. 

Remarks: Because of quantum backreaction effects, which are taken into ac- 
count implicitly, the deformation parameters are not necessarily constant during 
the evaporation process but in general some function (e.g. of the retarded time 
u). The first postulate implies that the action may change continuously with 
the retarded time u such that effectively one has to patch the classical solutions 

in dS and SS near r = 0. A minimal BH size of order of Planck mass is found and a cold 
remnant is reached after infinite time (note: a similar cold remnant has been found also in [14] 
for Einstein-Dilaton-Gauss-Bonnet gravity). 

^ Because of the rigidity result for consistent deformations of Einstein(- Yang-Mills) theory of 
gravity in _D > 2 [17] (at least) one of its premises must be violated if nontrivial modifications 
of the gauge symmetries are required to arise. In the present paper it is the dimension, D = 2, 
and the nonEinstcinian nature of the (dilaton deformed) nonlinear gauge symmetries. 
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of different models at lines u — const. It is not intended to construct one geo- 
metric model describing the BH evolution from the early stages until the final 
flicker. Rather, we have a continuous family of models at our disposal, labelled 
by the deformation parameters, and at each stage of the process one of them is 
the most adequate to describe geometry. 



Implementation: The crucial first condition restricts the possible deforma- 
tions of classical gravity due to quantum effects to consistent deformations. 
Classically, the SS BH can be derived from a Poisson-sigma model (PSM) with 
three target space coordinates [19], one of which can be identified with the 
so-called dilaton field (for the SS BH the dilaton field is proportional to the 
"surface area"). Of course, classical equivalence by no means implies quantum 
equivalence, but nevertheless the PSM formulation will be taken as starting 
point. Izawa has shown that the most general consistent deformation of a PSM 
is another PSM with the same number of target space coordinates [20]. Thus, 
the first assumption restricts to the class of PSMs with three target space coor- 
dinates. It contains infinitely many possible models, but with a little bit of extra 
structure each of them can be interpreted as a 2D dilaton gravity model [11]. 
For most discussions this first order formalism, which uses the Cartan variables 
as basic fields, is superior — for a comprehensive review cf. e.g. [22]. However, 
because more people are familiar with the second order formalism and since few 
technical details are needed the first order version of H2.1|l will not be employed, 
apart from the self-contained concise presentation in appendix^ The second 
order action reads [23] 



^SOG 



M2 



X§-^iS7Xr + ViX) 



(2.1) 



where X is the dilaton, g the 2D metric and R the Ricci scalar. The potentials 
U, V are arbitrary functions establishing the "[/ — V family of models" (cf. 
appendix^. It covers practically all relevant dilaton gravity models, including 
the CGHS [24], the Jackiw-Teitelboim model [25], the Katanaev-Volovich model 
[26] , and spherically reduced gravity [27] . By deforming the model in principle 
one can leave this class, but for simplicity it will be assumed that this does not 
happen. 



2.1.2 Second assumption: Asymptotic flatness 

Toy example: The "asymptotic" region hosting the observer is fixed by the 
asymptotic condition g(0) = 0. 



General form: The asymptotics of space-time has to be fixed in some way — 
typically by requiring asymptotic flatness during the whole evaporation process, 
but different behavior (such as asymptotic (A)dS) is conceivable. 



Specific form: The spacetime is assumed to be asymptotically flat during the 
whole process of evaporation ("stability of and J^"). 
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Remarks: This corresponds to a far field approximation — it is assumed that 
the observer measuring the Hawking flux is sufficiently far away from the BH 
and that this feature does not change during the evaporation process. It implies 
that at least for a distant observer the continuous patching of models is sensible, 
since all models have essentially the same asymptotics by construction. 

Implementation: Within the PSMs this assumption constrains us to the im- 
portant class of Minkowski ground state (MGS) theories, i.e. to those models 
which contain Minkowski space among their classical solutions. The MGS con- 
dition implies in the context of the a^b family defined below the linear relation 
a = 6 + 1. 

2.1.3 Third assumption: Avoidance of extremality 

Toy example: There are restrictions on the "causal structure" : the deformed 
potential in (|C.2|) is assumed to depend only on q and it must not have an 
extremum besides q = 0. As a technical simplification V (x q^ has been required, 
but supposedly this can be relaxed without changing the attractor solution. 

General form: Restriction of the causal structure (e.g. by requiring that the 
number of (non-extremal) horizons is non-increasing). 

Specific form: At most one (non-extremal) apparent horizon is present. For 
technical reasons only deformations within the so-called a — b family of dilaton 
gravity models will be considered [21], i.e. potentials of the form 

U{X) = -^, V{X) = -^X^+\ a,b,BeR, (2.2) 

to be inserted into the action H2.1(l . All of them exhibit at most one (non- 
extremal) horizon. 

Remarks: This is the most technical and from a general point of view the 
least essential assumption; it could be that it can be dropped in certain gen- 
eralizations, but for safety it is included here. It is needed to avoid extremal 
cases, both in the literal and the transferred sense. It can be considered as a 
consequence of the Penrose theorem (an isolated BH cannot split into two or 
more BHs, cf. e.g [28]), although of course one has to be carefuP by applying 
these classical theorems to BH evaporation (clearly Hawking's area theorem [30] 
is violated). 

Implementation: In this simplified scenario the only deformation parameters 
are a,i3 S M. For SS a = 1/2. The parameter B is almost irrelevant as it just 
defines the physical scale. Note that a, b appear nonlinearly while B enters only 
linearly. The generalization to the U — V family is treated in appendix IbI 

^E.g. in the Frolov-Vilkovisky model [29] although classically only one horizon exists at the 
quantum level an inner (apparent) horizon emerges. I am grateful to V. Prolov for enlightening 
discussions on this model. 
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2.1.4 Fourth Assumption: Boundedness of Hawking flux 

Toy example: The "toy model Hawking temperature" IjC.Sp remains constant 
during tlie wliolc dynamical process (actually this can be replaced by the much 
weaker condition of boundedness and positivity by a reparametrization of the 
asymptotic observer) . 

General form: The fourth assumption implicitly supposes that the notion of 
an asymptotic Hawking flux does make sense. While this is definitely true for 
the early stages of the evaporation it has to be regarded as a working hypothesis 
for the final evolution. Because this can be considered as controversial, before 
presenting the actual assumption some preemptive comments are in order: as 
mentioned above, we are not dealing with one model, but with a whole family. 
Each family member describes e.g. an eternal BH and as such allows for a 
sensible definition of Hawking flux and temperature. It is emphasized that at 
each instant the semi-classical approximation implicit in the notion of Hawking 
flux has to be fulfilled only for an infinitesimal amount of time. The crucial 
assumption is boundedness and positivity (or some stronger requirement) of 
the Hawking fiux as measured by an asymptotic observer. 

Specific form: An observer at is assumed to be isothermal during the 
whole evaporation process, i.e. the asymptotic Hawking flux does not change 
with the retarded time u (apart from the points where it is turned on and off). 

Remarks: Physically, this is the most important assumption because bound- 
edness typically implies that the ensuing dynamical system to be derived from 
these assumptions encodes effects beyond the semi-classical approximation and 
positivity ensures that no exotic matter flux is measured in the asymptotic 
region. The requirement of boundedness of the Hawking fiux is in the same 
spirit as the so-called Limiting Curvature Hypothesis (cf. [15] for an application 
to 2D BHs) which postulates that due to quantum gravity effects curvature 
remains bounded (thus implicitly assuming that the notion of curvature still 
makes sense). In that context one has to fix the potentials U,V in H2.1|) by 
hand appropriately. In this manner an eternally radiating remnant has been 
predicted [16]. In the rest of this paper the terms "asymptotic Hawking fiux" 
and "Hawking temperature" will be used interchangeably. 

Implementation: The fourth assumption is probably the least trivial one. 
In fact, it can be replaced by the much weaker condition that the Hawking 
temperature Th remains bounded and non- vanishing in a finite interval at . 
This is not unplausible from a physical point of view: it is supposed that quan- 
tum backreaction effects prevent the unbounded growth of the asymptotic fiux 
as predicted by semi-classical calculations; moreover, the Hawking temperature 
does not vanish up to the point where the evaporation terminates. The asymp- 
totic observer will be able to measure a temperature profile Th{u) as a function 
of the retarded time ii. Imposing that the Hawking process starts at ?i = = u 
by a regular diffeomorphism 




(2.3) 
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Figure 2.1: Example of a transformation to constant Th 



the new profile Th{u) is constant (see fig. 12.1(1 . It is convenient to rescale Th 
linearly such that it coincides with Th{u = 0) because in this way u and u are 
practically identical during the early stages of the evaporation process (dashed 
line in fig. 12.11) . Some physical interpretation relevant for the noninertial ob- 
server is presented in sect. 13.61 We recall that the evolution is understood to 
be within a continuous family of models, and not within a single model; thus, 
the diffeomorphism ((2.3|l is not a diffeomorphism acting on a single metric, but 
rather redefines the way in which the models are glued together. The knowledge 
of the "correct" profile Th{u) would be needed for determination of the precise 
form of the function f{u) in the Bondi mass, eq. (|3.5|) below. While the knowl- 
edge of f{u) would be necessary to answer questions like "how does the BH look 
like 10^ years after the evaporation starts?" the goal of the present analysis will 
be much more moderate: we are only able to determine the possible endstates 
of BH evolution. Thus, we are asking the following question instead: "With 
the given assumptions leading to a well-defined family of possible models and a 
well-defined evolution within this family, which are the attractor solutions?" As 
will be shown below it can be answered by exploiting isothermality Th = const. 



2.2 Exploiting isothermality 

For the a — h family of models obeying the MGS condition the Hawking tem- 
perature Th reads^ [31] 

Th = f{B,a){2MBH)'^''-^^'\ (2.4) 

where Mbh is the BH mass, a is the family parameter and f{B,a) is a pro- 
portionality factor. Th is related to the asymptotic Hawking flux Tg^^ via the 
Stefan-Boltzmann law in 2D: 

T^. = In . (2.5) 

At each instant the system behaves as if it were an eternal BH with a mass- 
to-temperature law as given by 1(2. 4|l in accordance with the discussion in sect. 
12.11 Infinitesimal changes of the temperature can be induced by changes of the 
black hole mass, changes of the family parameter (deformations) and changes 
of the normalization constant /: 

In i2MBH) + J (2.6) 



dTn 


cIMbh da 


a-1 


Th 


Mbh In (2Mbh) a(l - a) 


a 



'Natural units c = h = G m = kg = 1 will be used exclusively. Thus, A/pi^nck = 1- 
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It turns out that the first terms dominate over the Planck scale contribution pro- 
vided that Mbh is large enough. The term df / f is called "Planck scale contri- 
bution" because of the identity df = dTniMBH = 1/2) + 2f^dMBH{MBH = 
1/2). For large initial BH masses it can be neglected because the relevant terms 
in H2.6|l scale with In Mbh while the df-term is of order of unity. This is the 
reason why the scale parameter B is irrelevant until the final stage. For the 
moment this term will be dropped, but it will be reconsidered in in sect. 12.41 
In the rigid case df ^ = da the well-known relation (cf. e.g. eq. (6.29) of [22] 
and ref. [31]) 

dTn _ g - 1 dAlBH ,^ 
1h a Mbh 

is recovered (the choice a = 1/2 for the SS BH yields the famous inverse pro- 
portionality between mass and temperature). 

Isothermality dTn = implies (for a ^ 0, a 1, df = and Mbh 1/2) 
an ordinary differential equation which can be integrated from to Mf and 
from Qi to Uf, establishing 

_ln(2M/) ttij 
In (2M,) 1 - a^j 

which is already the main result of these considerations. 



2.3 The CGHS model as repellor 

Applying the previous analysis to a given starting model (e.g. = 1/2 for SS) 
with a given initial BH mass Mi ^ 1 one can predict with H2.8|l the "final" 
model at a given smaller BH mass. For z « 1 one obtains a/ ~ a^, as expected. 
Due to the logarithmic behavior of z the BH has to radiate most of its mass 
before the model is deformed appreciably (for instance z = 1/2 is induced by 
2Mf = ^/2Mi, which is just a tiny fraction of the original mass, but still large 
as compared to Afpianck)- The main questions is: what happens to af for 
monotonically decreasing z (especially in the limit z —> 0)7 It can be answered 
straightforwardly: the parameter af tends to zero unless > 1 (or = ±cx)). 
Note that z = corresponds to Mf — Afpianck/2. It will be assumed for the 
moment that this is the endpoint of the evaporation process. 

In language of dynamical systems the point a = is an attractor (a fixed 
point of the evolution which attracts all initial values which are close enough — 
in the present case close enough means Ui S (— cxd, 1)), while the CGHS model is 
a repellor^ (a fixed point of the evolution which repels all other initial values): 
each model with a = l— e, e>Ois driven towards the attractor, while each 
model with a=l-|-e, £>0 leads to a runaway solution a/ — *■ cxd because for 
small enough values of z a pole is reached in (|2.8|) . So ay = oo is an asymptotic 
fixed point which attracts all initial values € (l,cx)). The other asymptotic 
fixed point a/ = — oo is unstable. 

^It should be pointed out that for BH masses smaller than Mpianck/2 repellor and attractor 
change their roles: in this regime almost all models are driven towards the CGHS. But in this 
regime the contribution from df no longer is negligible (see next subsection). 
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2.4 Evolution close to the attractor 

Before addressing the further evolution of the attractor solution it is worthwhile 
to discuss it geometrically. At the 2D level the gauge symmetries of the model 
are related to a nonlinear (finite W) algebra and thus a geometric discussion 
in simpler terms would be desirable. To this end it is helpful to note that 
all models on the line a = b + 1 correspond to spherically reduced theories, 
where the dimension of the isometry sphere is given by 13 = 1/(1 — a) (so for 
a = 1/2 the 2-sphere corresponding to the isometry orbits of SS is obtained). 
In a sense, Hawking radiation induces an evaporation of all of the dimensions of 
the sphere but one, but the limiting solution in general does not correspond to 
"spherically" reduced Einstein gravity from 2+1 dimensions. Only for S = the 
model corresponds to a toric reduction of 2 + 1 dimensional Einstein gravity with 
the trivial line element (ds)^ — 2drdu + [duY — r'^idcj))'^. Note that it is not 
necessary to adopt the viewpoint of "evaporating dimensions" — alternatively, 
one can adopt a purely 2D point of view. It just seems easier to interpret 
the deformed gauge symmetries in terms of spherically reduced Einstein theory, 
because most people are accustomed to the gauge symmetries of the latter. 
The 2D line element (in Eddington-Finkelstein gauge) reads 

{dsf = 2drdu + (S In r - 2Mbh) {duf . (2.9) 

The curvature scalar R = —2B/r^ vanishes for r oo. The conformal diagram 
for H2.9|l is equivalent to the one for the SS BH [21]. 

For masses close to the Planck mass the a and B dependence of / cannot be 
neglected anymore and a refined analysis is needed. The more accurate version 
of H2.4|l reads 

T. = ^^(2M.,)(»-)/V (2.10) 

The scale parameter B has been fixed such that the next to leading order term 
of the Killing norm in an inverse "radius" expansion is normalized to —2Mbh- 
This seemingly harmless normalization has relevant consequences, so probably 
it is in order to say something about it: while there can be no debate about the 
leading order term (which can be chosen always equal to 1) the correct value of 
the next to leading order term is not completely trivial. The proposed normal- 
ization allows to identify Mbh at each instant with the mass of an eternal BH. 
This can be seen most easily by applying the standard subtraction procedure 
(cf. e.g. [31] or sect. 5.1 of [22]) 



M' 



= lim /?f(0 f - Vk{^) drX , (2.11) 



where K(r) is the Killing norm of the eternal BH with line element (ds)^ — 
2dudr + K(r)(du)^ and Kq is the Killing norm of the reference spacetimc (in 
the case of MGS models one takes naturally Minkowski spacetimc with Kq — 1). 
Applying (|2.11|l to the deformed SS BH (thus treating it at each instant quasi 
as an eternal BH) yields ikf'^s = Mbh- 

It is assumed henceforth that a lies in the relevant interval [0, 1) and thus 
X G [0, cx)). The refined variant of H2.8|l is given by 

^ , x,lnix/x,) Y' _ \n{2M) 

^ \n{2Mi) ) ' ln(2M,)- ^ ' 
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Figure 2.3: ...apart from the final flicker. 



Obviously, for large initial masses the correction terms are negligible (cf. fig. 
12.21 which displays the evolution of a SS BH with a moderate initial mass of 
1 gramme), as long as the ratio of xjxi does not become too small. However, 
close to the attractor this is precisely what happens. The BH evaporates down 
to Mpianck/2 already aX x — Xi (2Mi)-i/="', which is very close to the attractor. 

One can try to push the isothermality postulate to the very end. The BH 
eventually reaches a lower mass limit of 



2Miow = exp 



Xi , , , -l/xi 



e 



-(2Mi) 



(2.13) 



at Slow — (2Mj)~-'^/^'Xi/e. For a large SS BH this amounts to approximately half 
a Planck mass. From that point on isothermal evaporation actually increases 
the mass slightly until exactly Ajfpianck/2 is reached at the attractor solution 
(cf. fig. 12.31 which displays the final flicker of a SS BH with initial mass of 1 
gramme). Then only one question remains: is B in (|2.9|l equal or unequal to 
zero? A simple analysis proves that B scales with x close to the attractor, 
so the remnant geometry corresponds to (|2.9|) with B = 0. It is a regular 
geometry without any Killing horizon. Its "mass" just corresponds to the value 
of the Killing norm, which is globally normalizable to ±1 or 0, depending on 
its sign. In any case the geometry H2.9|l is flat for B ^ 0. For a positive sign 
of the last term the coordinates r, u = t — r have their usual meaning (radius 
and retarded time, resp.); for negative sign u corresponds to what usually is 
denoted by f = < + r (advanced time); for the vanishing case r becomes v. It 
is somewhat amusing that the mass has not really evaporated completely, but 
has just been converted into a purely geometric entity. 
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3 Physical discussion 
3.1 Approximate life time 

In this subsection (and only in this one) some actual meaning will be attached 
to the coordinate u. It will be assumed that the time it takes for the BH to 
evaporate as measured by the asymptotic observer is actually related to u. In the 
early stages this is obviously justified, because u ~ u, where u is the "physical" 
(retarded) time. However, for order of magnitude estimations the difference is 
negligible because as pointed out before the BH has to radiate away a lot of its 
mass before it is deformed appreciably. 

Applying the 2D Stefan-Boltzmann law H2.5|l to the energy loss of an evap- 
orating BH yields 

— cx - Th 2 = const. , 3.1 
at 

The proper life-time estimate then reads 

Uap (X (M,)"'+'/"' , (3.2) 

where the proportionality constant is close to unity. For a SS BH the approxi- 
mate evaporation time is given by tf^p « M?, which agrees with the standard 
estimation. The reason is that initially the BH for a long time looks like the SS 
BH — only in the last instances relatively dramatic changes do occur. 



3.2 The specific heat 

One of the extraordinary properties of the SS BH is its negative specific heat, 
being quite remarkable from a thermodynamical point of view. Thus it is of 
interest to check whether this feature changes or not. From (|2.10() we obtain 



dM 



BH 



dTi 



H 



l/(a-l) 

[^—ITh] . (3.3) 



1 - . xa/(a-11 / 1 — a 



Clearly, the specific heat is negative for the models considered here (a G (0, 1)). 
At a = 1/2 the SS result Cq cx is recovered. For a ^ the specific heat 
vanishes like a? /Th — thus, as the attractor solution is approached the specific 
heat vanishes. 

It has been advocated in refs. [32] that the negative specific heat is not 
only a sign of instability but also a sign of inconsistency of the canonical ap- 
proach. Applying a microcanonical analysis a mass loss rate different from the 
canonical one has been obtained. Thus, also the life-time estimate differs from 
(|3.2(l . It could be interesting to repeat the analysis of the present work for the 
microcanonical case. 



3.3 Semi-classical issues 

Since so many results are available at semi-classical level it is appropriate to 
discuss at least some of them from the perspective developed in this paper. 
First of all, the semi-classical derivation of the Hawking flux in the context of 
2D dilaton gravity (cf. e.g. [33]) and the corresponding mass-to-temperature law 
has been applied to obtain (|2.1U|) : in a sense, we have pretended that at each 
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step in the evolution the system can be described by an eternal BH at semi- 
classical level, but which of the infinitely many eternal BHs is the most proper 
one changes according to the evolution equation for the deformation parameter 
(|2.12() . So each of the intermediate "eternal" BHs has to fulfill the semi-classical 
approximation only for an infinitely small time interval. 

We should mention that spherically symmetric quantum deformations of 
the SS BH have been considered perturbatively already a decade ago by semi- 
classical means [34]. Also, the fact that dimensional reduction and renormaliza- 
tion need not commute does have a certain impact on semi-classical considera- 
tions [18]. Another interesting development is related to the question whether 
a near extremal (Reissner-Nordstrom, Kerr or deSitter) BH can evolve semi- 
classically towards its extremal limit. There is some evidence that the answer 
is negative [35] . This would imply that our restriction to non-extremal horizons 
is justified and singles out the attractor solution as the most likely endstate of 
BH evaporation in a more general frame. But clearly semi-classical considera- 
tions are not sufficient to decide these important questions conclusively. Due 
to the relation Th cx 1/Mbh semi-classical methods work very well only for 
macroscopic SS BHs. 

A complementary ansatz has been developed in [36,37]: instead of integrat- 
ing out the matter fields on a given background, geometry has been integrated 
out (exactly) and perturbation theory has been imposed upon the ensuing non- 
local non-polynomial effective theory, depending solely on external sources and 
matter degrees of freedom. This approach preserves unitarity and works very 
well only for microscopic (virtual) BHs [36,38]. Since the involved masses have 
to be much smaller than Mpianck the absence of remnants in these papers is not 
in contradiction with the results of the present work. 



3.4 The line element 

By a continuous patching of models along u = const, lines a global 2D line- 
element of an evaporating SS BH with initial mass Mi = Af adm can be presented 
as 



{dsf = 2drdu + {duf 



e{uf -u)[l 



j.x{il) 



+ 9{u~Uf)K^arn{r, u) 



(3.4) 



with x{u) as given by (|2.12|l (with z{u) = In (2MBondi(u))/ln (2M^£)Af)) and 

MBo„di(^^) = Madm - f{u)e{u) . (3.5) 

The Bondi mass is assumed to decrease monotonically (apart from the final 
flicker), i.e. f{u) is monotonically increasing and /(O) = for reasons of conti- 
nuity. The value f{uf) is determined by patching to K-ccra{'r, Uf) which encodes 
the last flicker. The specific form of f{u) depends on the model under considera- 
tion. The "remnant mass" Mrem is of order of Mpianck, but the further evolution 
depends on K^caiir^u); it has been argued in sec. 12.41 that it will look like the 
Killing norm in 1)2.9(1 with B = Q. The explicit solution for x{u) involves the 
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Figure 3.1: Apparent horizon as a function of AfBondi 



Lambert W function® {W{z) is a multivalued function exhibiting algebraic and 
logarithmic singularities, but above the limit H2.13|l only the principal solution 
for W{z) in z — W is needed — for masses below A/pianck/2 two real branches 
exist, but a simple continuity argument establishes one of them as the proper 
one; thus, no branch cut ambiguities arise, except for the point where the limit 
(I2.13|) is saturated): 

2Madmx{u) = \\ , m{u) := 2Madm In (2A/Bondi(M)) • (3.6) 

For a Bondi mass below Mpianck/2 a second real branch of W exists; however, 
continuity of the evolution restricts us to the principal branch. The apparent 
horizon is located at r{u) = rh{u) with 

Thiu) ^exp[W[m{u)]]. (3.7) 

It is a monotonically decreasing function with u (see fig. I3.1|l . If we had 
taken the simpler eq. (|2.8|) instead of (|2.12|) the horizon were located always 
at fh = 2Madm — hence, the Planck scale contribution df in H2.5|l is responsible 
for a decreasing apparent horizon. At the lower mass limit H2.13() the horizon is 
located at r™™ — 1/e. From that moment on, as pointed out before, isothermal- 
ity would imply an increase of A/Bondi, somewhat in contradiction to common 
sense. 

The corresponding Carter-Penrose diagram is depicted in a sketchy manner 
in fig. 13.21 It has to be taken cum grano salis. For simplicity, region I is taken 
as fourdimensional Minkowski spacetime. At v — vq a light-like shock wave 
forms the SS BH. Region II is the SS region (in this region eventual radiation 
from the formation process is emitted to which has been blithely ignored). 
At M = Hawking radiation commences (the light-like line at m = would 
be the SS horizon if no evaporation occurred). In region III the continuous 
deformation towards the attractor solution happens and the apparent horizon 
decreases by virtue of 1)3.7(1 . At u = Uf the minimal mass A/iow is reached. The 
last flicker happens approximately at u = m/ (this line has in fact a width of 
order of Planck length; thus also the full dot at the origin is in fact not a point 
but a region of order of Planck area). Region IV is the remnant geometry, i.e. 
according to the previous discussion it is Minkowski spacetime. 

^For a comprehensive discussion of W{z) see [39]. I am grateful to A. Rebhan for providing 
this reference. 
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Figure 3.2: A sketchy conformal diagram depicting the evaporation of a SS BH 



It is instructive to calculate the 2D curvature scalar 

2a(u)MBondi(w) 



R{r, u) 



(1 - a(u))2r(2-'^("))/(i-"(")) 



(3.8) 



For SS the result R = —4M/r^ is recovered. Of course, the fourdimensional 
curvature scalar vanishes for SS, but not all scalar invariants built from the 
Weyl tensor do. They are encoded in (|3.8|l . 



3.5 Speculations on the last flicker 

There are three simple ways to deal with the very last stage: either suppose that 
a small amount of exotic matter really does exist; then the attractor solution is 
the endpoint of evaporation. However, since one of our postulates was actually 
positivity of the asymptotic flux this violation, as moderate as it might be, is 
not a pleasant feature. 

Alternatively, it seems tempting to give up isothcrmality as the lowest mass 
is reached to avoid exotic matter, i.e. to take the line element (|3. 413.6(1 . but 
assume that (in spite of fig. 12. 3|) the Bondi mass is decreasing below Miow as 
given by (|2.13() . This has some drastic consequences: x(u) suddenly acquires 
a non-vanishing imaginary part, x{u) = a + i(3, with a,/3 G R. Thus, the line 
element stays real only at the discrete set of points r = with n e Z and 
ro = exp [tt//?] (note the duality vq — s- l/ro). There is a Z2 grading depending 
on whether n is even or odd (the quantity r^^ becomes positive or negative, 
respectively). Eventually the real part goes to zero: 

2AfBondi(Mo) - exp [-Tr/{U'lADM)] < 2M10W x{uo) = i(3 . (3.9) 

The quantity /3 is positive and in fact very small close to the critical point. 
One can only speculate about the interpretation of the appearance of a dis- 
crete structure (or alternatively the appearance of a complex metric). On a 
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superficial level complex and/or discrete structures are expected to arise in a 
comprehensive theory of quantum gravity; maybe it is most proper to interpret 
this effect as a phase transition: above Mjow an adequate description can be 
achieved by a (smooth) metric, but at il/iow the system "condensates" to some 
discrete/complex entity. It could be worthwhile to investigate this idea more 
rigorously^. Applying the results for the asymptotic limit of the Lambert W 
function, W{z — )■ oo) = In |z| + . . . where the omitted terms contain multiple 
logarithms and branch cut ambiguities [39] , allows a discussion of geometry in 
the hmit 2AfBondi =: £ ^ 0: 

limfl-— — —j] = 1- lime!-''^'-/''^!'"^! =1 (3.10) 

This produces just the attractor solution. 

The third possibility (which is the one favored by the author because of its 
simplicity and because of the issues discussed below in sect . 13.611 is to assume that 
the evaporation process simply stops as the mass drops down to Mpianck/2. This 
corresponds to a discontinuous jump from the principal branch to the second 
real branch of the Lambert W function. Thus, neither a phase transition to a 
discrete structure nor exotic matter is needed. The "undershooting" as depicted 
in fig. 12. 31 has been eliminated by cutting the evolution process below z = and 
connecting the two (otherwise disconnected) points x{z = 0). One can call such 
a procedure also "phase transition" because the Hawking temperature jumps 
from a finite value to zero. 

Thus, all three approaches yield the attractor solution (|2.9|l with B = 
as endpoint — they differ "only" in interpretational issues: either one has to 
allow for (a tiny amount of) exotic matter or for a discrete/complex space- 
time structure or for a discontinuity as Afpianck/2 is reached. The second and 
the third possibility imply the occurrence of a phase transition, because both 
induce a non-smooth change at a branch point. The order of this transition 
differs: while for the second scenario continuity has been maintained (and thus 
one can call it a second order phase transition), the last version implies a first 
order transition. 

It has been shown in sect. 12.41 that the mass is effectively converted to a 
purely geometric quantity of the attractor geometry, but energy conservation 
dictates that something has to happen with the final Mpianck/2. There are two 
possibilities: 

1. A remnant in the form of matter concentrated at r = exists. 

2. A thunderbolt emits a total energy of Mpianck/2 at the end of the evapo- 
ration process. 

Of course, it could be that both possibilities have to be combined. 

Suppose for a moment the matter remnant scenario. A point particle in 
2-1-1 dimensions corresponds to a cosmic string in 3-1-1 dimensions. Then the 
question arises: Whatever happened to spherical symmetry — how can it be 
that a spherically symmetric object (the SS BH) evolves into an object with a 

^ Since we are already speculating it should be mentioned that the effect of conifold sin- 
gularities in string theory [40, 41] appears to be similar (even the function x\nx plays an 
essential role in both contexts, cf. 12.121 in the current paper and (2.4) in [40]). I am grateful 
to M. Kreuzer for discussions on that topic. 
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preferred axis (the cosmic string)? It should be emphasized that the total mass 
of the cosmic string is bounded by Afpianck/2; thus, if it were infinitely long its 
mass density per unit length went to zero. Possible interactions between cosmic 
strings and SS BHs have been discussed in the context of "thorny spheres" in 
ref. [42]; for a physical discussion see section 7.4 of [1] and references therein. In a 
more physical setup including also angular momentum it is conceivable that the 
axis of the remnant string is aligned with the axis of rotation of the macroscopic 
BH. Moreover, its length is unlikely to be infinite — a typical length scale to be 
expected would be around Planck length and thus also the mass density will 
be of order of unity in Planck units. Such a finite size object violates spherical 
symmetry only moderately. 

3.6 Dilaton evaporation 

As mentioned the variant of a first order phase transition is favored by the 
author. While revising this paper a new work appeared [43] which supports 
it. When talking about a phase transition one would of course like to know 
which are the two phases and what separates them — after all, so far we have 
achieved only a smooth transition towards the attractor geometry, so what can 
possibly provide the degree of nonsmoothness pivotal to a first order transition? 
The answer is provided by the dilaton field: while before the evaporation is 
terminated the dilaton necessarily is non-constant (it can be rearranged to be 
linear in one of the coordinates by exploiting the gauge degrees of freedom), 
afterwards it can be constant. Therefore, the phase transition could be from 
a linear dilaton vacuum to a constant dilaton vacuum. It will be argued in 
this section that indeed this is the case. Thus, the current approach provides a 
realization of the semi-classically motivated conjecture [44] that constant dilaton 
vacua might be the final state of BH evaporation. 

Most dilaton theories H2.1|l allow for one or more constant dilaton vacua 
X = const., where the constant is determined by the zero(s) of the potential® 
V{X) (cf. sect. 2.1 of [45]). The corresponding geometry can only be (A)dS, 
Rindler or Minkowski space. Thus, the only constant dilaton vacuum exhibiting 
the MGS property is Minkowski space. Because also the attractor solution (|2.9() 
with i? = is Minkowski space it seems that the dilaton has evaporated to a 
constant dilaton vacuum. However, there is an important subtlety (cf. sect. 3 
of [45] ) : if the attractor solution is patched to such a constant dilaton vacuum 
along a u = const, line then a matter flux is induced on it, which is remarkable 
insofar as curvature is continuous (namely zero in our particular case). This 
provides a mechanism for a shock wave like evaporation of the final Mpianck/2. 

At this point it should be recalled that the observer is a noninertial one — 
naive semiclassical extrapolation implies that the boost velocity has to approach 
the vacuum velocity of light in the limit of vanishing mass; this is nothing but 
an ultrarelativistic (Aichelburg-Sexl [46]) limit! Incidentally, the shock-wave 
induced by such a limit in 2D dilaton gravity takes the same form as the shock- 
wave which arises when patching a constant dilaton vacuum to a non-constant 
dilaton vacuum [43]. Thus, the scenario with the first order phase transition 
fits nicely together with the boost of the observer implicit in the assumption of 
isothermality. Actually, exploiting the results of [43] one can even provide an 

®Note that for the attractor solution V = and thus the dilaton is unconstrained. 
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educated guess for the total energy emitted in the thunderbolt: only half of the 
final BH mass will be emitted, because the discontinuous directional derivative 
of the dilaton, in the notation of [43], jumps from a constant to minus 
this constant. However, patching to a constant dilaton vacuum implies that it 
jumps only to zero, i.e., half of the amount as compared to an ultrarelativistic 
boost — because the height of the jump sets the scale for the energy emitted by 
the shock wave it seems that only half of the mass of the final state can be 
emitted in the thunderbolt, the other half will reside in a cold remnant. 

Although the first order formulation so far has been avoided it will be con- 
sidered in this paragraph because there is a nice target space interpretation 
of the phase transition. A reader who is not familiar with PSMs can skip this 
paragraph or may consult appendix^and references therein. There are two ba- 
sic differences between symplectic manifolds and Poisson manifolds: the Poisson 
tensor can have a nonvanishing kernel and its rank need not be constant. Indeed, 
in dilaton gravity there is always a nonvanishing kernel and the corresponding 
conserved quantity is related to the total energy (essentially the ADM mass, 
whenever this notion makes sense). But the rank is constant, namely 2, except 
for eventual isolated points where it vanishes (like the bifurcation 2-sphere in 
the SS BH). The evaporation to a constant dilaton vacuum corresponds to a 
transition to a region where the rank of the Poisson tensor vanishes! Thus, 
before the phase transition the world sheet geometry is defined by symplectic 
cuts in target space, while afterwards it is defined by a single point in target 
space. 

To conclude, the phase transition is not one between two different geometries 
(because the line element for the attractor solution is equivalent to the one for 
the constant dilaton vacuum), but between different behaviors of the dilaton 
field. It can be understood in terms of an ultrarelativistic boost of the observer. 
A similar mechanism could be of relevance for scalar-tensor theories in 13 = 
4; there, a constant dilaton vacuum simply corresponds to a phase where the 
effective Newton constant does not vary. 

3.7 The information paradox 

It is difhcult to avoid addressing the information loss puzzle (cf. e.g. [47]) when 
discussing the long time evolution of BHs. The main two alternatives are 1. the 
information is lost, as proposed by Hawking [48] and 2. some mechanism exists 
that prevents the information loss. Most people with a field (or string) theory 
background favor the second variant. In the following it will be discussed very 
briefly which of the proposed solutions fits into the current approach. 

The quantum hair solution does not fit into the current formalism since the 
formation process has been neglected. Also the possibility that the information 
is emitted together with the Hawking radiation is not applicable here, at least 
not without specifying the matter content of the theory; neither is the baby 
universe scenario, at least not in an obvious way. Since a remnant geometry 
has been predicted there remain two natural solutions: either the information is 
retained in the remnant or it emerges at the end of the evaporation process. Due 
to the fact that the proposed remnant geometry is flat it is hard to see where 
any information could be stored (apart from the scenario with a cosmic string 
remnant), so the only possibility seems to be that the information is released 
during the final stages. Indeed, this is what the Carter-Penrose diagram fig. 13.21 
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seems to suggest, as well as the discussion in sect. 13.61 It will be studied in more 
detail below. 

The information paradox manifests itself typically in a pathological feature of 
Carter-Penrose diagrams related to BH evaporation (cf. e.g. fig. 1 of [49] for the 
classical example; in that reference it is argued in favor of a vioation of locality). 
Thus, it is worthwhile to study the present Carter-Penrose diagram fig. 13. 21 and 
to check what, if any, pathological features are encoded there. To this end the 
consideration of test particles will be useful. Note that in the language of the 
S-matrix approach to quantum black holes [50] the "hard matter" components 
implicitly have been taken into account by our assumptions. "Hard" refers to 
the fact that its backreaction upon geometry cannot be neglected, as opposed 
to "soft" matter which can be treated like test particles. It will be assumed for 
simplicity that the test particles are massless. 

Let us start at some value of advanced time v = Din at and follow the 
light rays until they are reflected at the origin and eventually emitted to a point 
at corresponding to a certain value of the retarded time u = Uout- If ^^in 
is sufficiently close to i~ then Uout < 0, i.e. the outgoing ray will pass solely 
through region I, encounter the infalling matter flux and finally pass through II. 
Similary, for values of Wjn sufficiently close to i'^ the lightray will pass through 
regions II and III, encounter the outgoing shock wave and end up in region IV 
where eventually it will be emitted to J^'^ close to Actually, this is true 
for all values > vq. Of course, rays starting closer to vq will additionally 
pass the apparent horizon, but so far no pathological features were found. At 
V = Vq there is the problem of the potential singularity encountered at the 
origin. Recently it has been proposed to impose a boundary condition at the 
singularity [51] (cf. [52] for a critical comment) as a possible path to a solution 
to information loss; in the present context it is not quite clear whether the 
singularity exists at all as it lies in the very region we could only speculate 
about. So let us diregard it and focus on the region at that has not been 
covered so far: Vm < vq but still sufficiently far away from such that light 
rays are not scattered to region II. 

Such a test particle will trespass the line u — and be reflected at the origin 
still in region I. However, it crosses the infalling matter shock wave between 
u = and u = Uf thus entering region HI. It will follow on a curved path more 
or less parallel to the apparent horizon and then pass the outgoing shock wave. 
In region IV it moves first again towards the origin and after the reflection it 
flnally reaches J^+. Such test particles are remarkable for three reasons: 1. they 
are reflected twice at the origin, 2. they cross infalling and outgoing shock wave, 
3. for each test particle travelling in this way there exists exactly one "partner 
particle" emitted at Win > vo- The "partner particle" will end up at the same 
point at J^"*". Thus, tracing back the light rays from J^"*" there exists a region 
between u — Uf and some larger (but finite) value of u where the light ray splits 
into two as soon as the outgoing shock wave is encountered from above: one 
directly going back to ^~ and the other one entering region HI. In this sense, 
the outgoing shock wave causes birefringence. The line w = acts like a mirror 
regarding the two "partner particles" : if one of them is close to (far from) vq so 
is the other. 

On top of that, the later the test particle is emitted at in region I (i.e. 
the closer Uin is to vq) the earlier it arrives at J^"*" in region IV (i.e. the closer it 
is to u = u/). Thus, causality is violated (or rather: inverted), but this violation 
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is not noticed at ^+ until the BH has evaporated completely! 

This is the pathological feature we were looking for (or hoping to avoid): 
there is something very strange about J^"*": while the regions u < and close 
to are unproblematic, there is no way to send information to the region 
M/ > u > by means of test particles — this coincides with the region dominated 
by Hawking radiation caused by "hard matter" . Moreover, there exists a region 
Ue > u > Uf with some finite Ug where information may arrive in two ways: 
either directly from with only one reflection or with two reflections and 
inverted causality, as described in the previous paragraphs. It is worthwhile to 
study these features in more detail, but such a study would lead beyond the 
scope of this paper. 

4 Relaxing the conditions and conclusions 

RelELxations within 2D dilaton gravity There are three generic possibili- 
ties regarding the flrst assumption: first, one can assume the complete absence 
of deformations with the usual ensuing problems; second, one can allow for con- 
sistent deformations only as in the present approach (with the important caveat 
that only deformations of the spherically reduced theory have been considered); 
third, one can allow for arbitrary deformations which enlarges the "phase space" 
of possible models but at the same time reduces the predictability. It is sound 
to start with the second possibility until some fundamental theory dictates the 
precise form of the allowed deformations. 

Relaxing the second assumption is possible in principle — e.g. a BH immersed 
in (A)dS or Rindler spacetime can be considered along these lines. However, the 
fact that asymptotics is not changed essentially during the evaporation process 
should be maintained, because otherwise the notion of a long-time asymptotic 
observer will be challenged. For the third assumption similar considerations 
apply: one can allow e.g. for 2 Killing horizons and thus is able to discuss the 
Reissner-Nordstrom case. Such a generalization is not as harmless as it might 
seem at first glance: once 2 or more horizons are allowed extremal ones may 
emerge during the evaporation process. Since they do not produce Hawking 
radiation extremal geometries are additional candidates for possible endpoints 
of the evolution. Because stringy methods are capable to describe (nearly) 
extremal BHs very well (cf. e.g. [53]) one could imagine that the non-extremal 
evolution is described as discussed in the current work and from the extremal 
point on string theory takes over. However, there is classical and semi-classical 
evidence that nearly extremal BHs do not evolve towards the extremal limit, 
which can be taken as a justification of our neglection of extremal horizons (cf. 
sect. 13. 3|) . 

Relaxing the crucial fourth assumption does not seem plausible from a phys- 
ical point of view — either the asymptotic flux would have to be negative at some 
point ("anti-evaporation") or it would have to grow without bound. 

Obviously it is desirable to generalize from the a — b family to a more general 
class (or the most general class) of possible models. Some steps in that direction 
are presented in appendix^ As long as no extremal horizon is present evidence 
points again towards the previously discussed remnant geometry as the endstate 
of evaporation. 
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Beyond 2D dilaton gravity The restriction to PSMs^, even though it allows 
for phenomenologically relevant models like the SS BH or generalizations thereof 
(with charge and/or cosmological constant), could be overcome by leaving the 
comfortable realm of 2D. Due to the rigidity result [17] mentioned in footnote 
13 a natural framework would be a scalar-tensor theory a la Jordan-Brans-Dicke 
[55] or its more recent incarnation as quintessence [56]. A brief discussion of 
scalar-tensor theories in the framework of 2D models can be found in [57]. 
It should be possible to apply the methods of this work to the dimensionally 
reduced models treated in that reference. 

To be a bit more concrete: an action of a I?-dimensional scalar tensor theory, 



contains two arbitrary functions U, V of the dilaton X which usually are ad- 
justed either by hand or by some guiding principle from a more fundamental 
theory. Tinkering with the assumptions, it should be possible to generalize the 
present work from D — 2 to arbitrary D, in particular to D = 4. If again 
the phenomenon of "dilaton evaporation" to a constant dilaton vacuum occurs, 
the final state will be a solution of ordinary vacuum Einstein theory with a 
cosmological constant A, the value of which depends on the function V{X) of 
the attractor solution: A oc V{Xcdv)/Xcdv, where Xcdv is a solution of 
the equation 2V{Xcdv) = XcdvV [Xcdv)- The existence of such a solution 
poses nontrivial constraints on the dilaton potential V . 

What are the experimental consequences? Since for each SS BH (in- 
dependent of its original mass, as long as it was larger than Mpianck/2) a 
cosmic string with a total mass of maximal Afpianck/2 remains this puts an 
upper bound on the total amount of primordial BHs. Actually, this bound is 
well-known [58] because several models predict a relic with a mass of order of 
A^pianck [59]. Apart from their total mass cosmic strings contribute to density 
fluctuations which should be visible in the observed cosmic microwave back- 
ground anisotropics (cf. e.g. [60] for an analysis of the BOOMERanG data). 
Finally, the thunderbolt with a total energy of maximal Mpianck/2 should lead 
to observable consequences, possibly in the form of 7-ray bursts. The sum of the 
energy of the thunderbolt and the mass of the cosmic string remnant is exactly 
A^pianck/2 provided energy conservation holds. However, none of these conse- 
quences provides a "smoking gun signature" . Moreover, as mentioned before, it 
is not unlikely that eventually non-spherical modes will become relevant during 
the evaporation process, which further limits the phenomenological relevance of 
the present work. 

Brief summary Let us finally repeat the lines of reasoning underpinning the 
approach presented in this paper: it is generally expected that during the final 
stages of BH evaporation some drastic changes will occur, since the semiclassical 
theory extrapolated to the very end predicts its own failure. In order to discuss 
the issue of long time BH evaporation in more detail either quantum gravity 
must be applied or some trick must be used. It has been proposed to follow 

generalization to graded PSMs in order to describe dilaton supergravity models should 
be straightforward, cf. [54], 




(4.1) 
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the latter route by assuming an isothermal observer at i'^, i.e. the asymptotic 
Hawking flux (and hence the Hawking temperature) as measured by this ob- 
server does not change during the whole evaporation process. It can be replaced 
by a much weaker condition, namely boundedness and positivity of the asymp- 
totic Hawking flux. As a consequence — together with "natural" assumptions 
regarding the asymptotics, the causal structure, the nature of allowed quan- 
tum deformations of gravitational symmetries and some simplifying technical 
assumption further restricting the allowed deformations — an evolution of the 
SS BH towards a remnant geometry (|2.9() with B = has been predicted. We 
speculated about the final flicker, argued in favor of a first order phase transition 
to a constant dilaton vacuum, discussed briefly implications for the information 
paradox and concluded with possible generalizations of our postulates. 

If nothing else, the present approach shows that quantum induced deforma- 
tions of the gauge symmetries of gravity can play a pivotal role for the under- 
standing of the long time behavior of black holes. 
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A Deformations of dilaton gravity 

This appendix is meant to be self-contained. Its purpose is to review concisely 
some features of ordinary and deformed dilaton gravity which are (at least 
implicitly) relevant for the rest of this paper. 

To this end (only in this appendix) the first order version of 12.1|l [19,61] 

i'^^ = / [Xa{D Ae f +Xd/\uj + eV{XaX'',X)] , (A.l) 

will be employed where X is the dilaton field, e° is the zweibein one-form, e is the 
volume two-form. The one-form lu represents the spin-connection w^f, = e^bW 
with Eab being the totally antisymmetric Levi-Civita symbol. The action l)A.l|l 
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depends on two auxiliary fields X°-. It is a special case of a PSM [19,62] 



L 



PSM 



2 



(A.2) 



IM2 

with the Poisson-tensor (/ and C are arbitrary functions of A"') 

which fulfills the (generalized) Jacobi identity 

P^P^'^ +cycL(z, J, fc) = 0, (A.4) 

and a three dimensional target space the coordinates of which are A, A*^. In 

light-cone coordinates (ry^ = 1 ~ ?/ |_, r;4--|_ = = f] ) the first ("torsion") 

term of l|A.Hl is given by 

XaiDAe)" = riabX\D Aef = X+ {d - to) A er + X- (d + lu) A e+ . (A.5) 

The function V is an arbitrary potential depending solely on Lorentz invari- 
ant combinations of the target space coordinates, namely A and A+A^. To 
match the first order with the second order formulation one has to choose 
V = ;7(A)A+A- + V{X) with the same functions U, V as in (EH. To match 
the PSM with the first order version one has to choose / = 1//(A), C = 
/(A)A+A- +w{X) with 

/X pX 
U{y)dy, w{X):=j I{y)V{y)dy . (A.6) 

Variation of the fields yields the equations of motion 

dX' + P'^Aj = , (A.7) 
dAA.-^P/'AkAA, ^0. (A.8) 

Under the symmetry variation (note that e = £(2:^, A*), where x'^ are the world- 
sheet coordinates and A* the target space coordinates) 

SeX' = P'%, , S,Aj = -de, - P/-"e„Afe , (A.9) 

the action (|A.2|I transforms into a total divergence 



S,L= / dX'si. (A. 10) 

JdM2 

The commutator of two symmetry variations closes only on-shell in general: 

[6,„S,,]X'^S,,X\ (A. 11) 

[S,,,Se2] = 6,,A, + {dX^ + P^'^Ak) P,."™ei™e2n , (A.12) 

with 

e3^ = P/'=eifc£2, + P^' [efc, (e.OJi^ , (A. 13) 
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and [/(£),5r(e)]i2 := /(ei)g(£2) - /(e2)5(£i)- 

The introduction of the Schouten-Nijenhuis bracket [63] 

{X\X^} = {^',»} = P'^dji*) , (A. 14) 

allows to reexpress e.g. the conservation law as {X\C} = 0, V«. According to 
Izawa [20], consistent deformations in the sense of Barnich and Henncaux are 
those which leave the bracket structure ljA.14(l essentially intact, i.e. the Poisson 
tensor may change its functional dependence on the fields X', but its dimension 
(and the fact that one stays within the realm of PSMs) remains fixed. 

Since Izawa's result is of some importance for the present approach his proof 
will be sketched briefly: starting point is 2D BF theory (the conventions of [20] 
have been translated accordingly) 

So= [ dX'AA,, (A.15) 

Jm2 

which is just a very trivial special case of a PSM ljA.2|l with vanishing Poisson 
tensor. Then, the minimal solution S to the classical master equation and 
the generator of BRST symmetry s are constructed. A deformation Li of the 
Lagrangian must obey [5] 

sLi+dai^^O. (A.16) 

This, together with its descent equations sai + dao ~ 0, sgq = 0, allows to 
construct the most general consistent Li. It is given by eq. (14) of [20]. The 
deformed action 

S = S+I Li=Sa + \( P'^ Aj A A, + antifieldterms, (A.17) 
Jm2 2 Jm2 

satisfies the classical master equation and reduces to (jA.2|) for vanishing anti- 
fields. Actually, it even satisfies the quantum master equation as there is no 
contribution from renormalization [64]. No further deformations occur because 
P'-* fulfills the (generalized) Jacobi-identity (|A.4|I . i.e. the first order deforma- 
tion of abelian BF theory, which is nothing but a PSM, provides already the 
most general consistent deformation. Thus, also the most general consistent 
deformation of a PSM is just another PSM (with the same number of target 
space coordinates). 

It is important to notice that dilaton gravity, if formulated as a PSM, requires 
additional structure, namely a correspondence between the gauge fields and 
the line element. The suggestive notation Ai = (w, e~, e"*") together with g = 

277-1 e+e~ apparently fixes this structure in the undeformed case, but there 

is no a priori relation between the three gauge fields and the line element if 
deformed PSMs are considered. This subtle and important issue is addressed in 
detail in [11]. 

B The U -V family 

For technical reasons the main part of the paper concentrated on the a — b 
family. In this appendix the results will be generalized to the U —V family, i.e. 
to models with a Lagrangian 1)2.1(1 with arbitrary potentials U{X) and V{X) 
which fulfill the MGS condition and which have at most one horizon. 
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All local and global classical solutions of such models have been discussed 
extensively in a series of papers by T. Klosch and T. Strobl [65]. For a first 
orientation on classical, semi-classical and quantum dilaton gravity in 2D the 
review ref. [22] may be consulted. 

In the following the notation and results of [37] will be used; for sake of self- 
containment we recall two definitions: the "integrating factor" I{X) and the 
conformally invariant combination w{X) as defined in l|A.6|l . The MGS property 
reads w oc 1//. Provided that the equation I ~ c has exactly one solution for 
positive values of c and no solution for c < the existence of at most one horizon 
is guaranteed. The constant c scales like 1/Mbh- The Hawking temperature 
reads Th oc U{Xh)MBH where Xh is the solution of I{Xh) = 1/Mbh- It is 
convenient to introduce the inverse function of the integrating factor, I* ol ^ id, 
because X^ — I* {1/Mbh)- With these preliminaries the Hawking temperature 
for MGS models of the U ~V family with one horizon reads 

I \i/Mbh 

Now we would like to address the issue of boundedness of Th ■ Suppose that 
/* is a given function of X (from this one can deduce all other quantities). First, 
it will be assumed that /* has a Laurent series expansion of the form 

N 

r ^ a„X", N el. (B.2) 



If iV = then the leading term drops out after differentiation, so in this case 
the value for N relevant in all subsequent formulas is given by the next to 
leading order term. For small values of the BH mass the Hawking temperature 
is dominated by the highest power N in l|B.2p (and hence the discussion is 
essentially reduced to the a — b family) , 

Th{Mbh ^^)~a{MBHf+\ aeM+. (B.3) 

Thus, only for N > —1 (but N 0, because this term drops out after dif- 
ferentiation) the temperature remains bounded for all times. Relevant exam- 
ples/counter examples are the CGHS model with iV = — 1 and spherically re- 
duced gravity from D — 4 with N — —2, respectively. The attractor solution 
does not fit into this scheme because / = const, does not allow for an inverse 
(but this model is rather trivial anyway). 

If there is a nonperturbative prefactor in front of the sum in (|B.2|I which 
is regular at Mbh = then again the highest term dominates and the same 
bound > — 1 is recovered. If there is a different (but at Mbh = regular) 
nonperturbative prefactor for each term then still the same conclusion holds. 
Thus, the bound A^ > — 1 is rather general. If the prefactor is singular a 
different bound may be obtained (e.g. A^ > — 1 for a logarithmic prefactor). 

If the sum in (|B.2|) does not terminate (i.e. A^ +oo) then an essential 
singularity at Mbh = is encountered. In order to decide whether Th remains 
bounded again the limit Mbh +0 has to be considered. In cases where the 
Laurent series can be resumnicd this is straightforward (e.g. for something like 
exp [-1/Mbh])- 
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In principle one can also extract an MGS theory for a given temperature 
law (if, for instance, the temperature law can be motivated by other means). 
The only MGS model with Th = const, turns out to be the CGHS. The only 
MGS model with Th = is the (trivial) attractor solution. Thus, if Hawking 
radiation is assumed to stop at a certain point the attractor solution is the 
unambiguous remnant (unless, of course, one allows for extremal horizons). 

A generalization of these considerations to non U — V theories (like e.g. the 
class of models discussed in [66]) is possible in principle, but somewhat tedious 
in detail. 

C A toy model from classical mechanics 

Suppose a dynamical system which can be described by a Lagrangian action 
depending on two fields, one of them being called "geometrical" (q) and the 
other one "matter" (Q). Suppose for simplicity that we are in 1+1 dimensions 
and that the geometrical part alone can be cast into a 0+1 dimensional form 
(assuming staticity), but that we do not know how to solve the coupled system. 
Instead, deformations of the geometrical subsystem will be studied and assump- 
tions will be made which finally allow to obtain an attractor solution without 
having to specify the matter part. 

To keep it basic, let the undeformed geometrical action be a harmonic oscil- 
lator 



where q :— dq/dr and r is a quantity to be interpreted as "radius". There is 
no gauge degree of freedom and one physical degree of freedom. In analogy to 
appendix^we want to find all consistent deformations in the technical sense of 
Barnich and Henneaux [5] — i.e. the number of gauge degrees of freedom is not 
changed as well as the number of physical degrees of freedom, but the action 
may change as well as the gauge symmetries. Often a perturbative analysis in 
terms of the deformation parameter is employed, but in the present paper the 
model is simple enough to go beyond perturbation theory. For the toy example 
deformations are rather trivial,^" 



where e is the deformation parameter and V{q, q, g, . . . ; 0) = 0. We will sim- 
plify the issue further by requiring independence of q or higher derivatives and 
impose linearity in q^ and in the deformation parameters (these requirements 
correspond to some of the "physical" assumptions that are imposed on the sys- 
tem and of course their validity should be justified; however, as this is merely a 
toy example this is omitted). Thus, we have a 1-parameter family of deformed 
geometrical actions at our disposal. 



^'^First, an antifield q* is defined with antibracket {q,q*) = 1; then the BRST charge is 
constructed, s = ijS/Sq*; finally, the deformation equation sL'^' + dag = together with its 
decent equation sao = have to be solved, where L^^^ is the deformed Lagrangian. These 
methods [5] (cf. also [6] and references therein) are being used in appendix E] but appear to 
be a slight overkill for the toy example. 




(C.l) 




(C.2) 



(C.3) 



25 



REFERENCES 



Initially, s — 1. Let us impose the "asymptotic" condition q{0) = and assume 
that the observer is located in the "asymptotic region" r = 0. Then, the general 
solution of the equations of motion to be derived from ljC.3|l is 

g(r;e ^ 0) = Asin(Ver-) , q{r; e = 0) ^ Br , A,BeR. (C.4) 

The asymptotic observer patches together these solutions along t = const, lines, 
where t is to be interpreted as "time" . Now let me define a "physical observable" 
which I would like to call "toy model Hawking temperature" : 

Th :=g(0) = Ax/i. (C.5) 

Assume that the system can be reparametrized in such a way, that in the full 
dynamical system there is no time-dependence of Th- Provided A is positive and 
strictly monotonically increasing, e has to be strictly monotonically decreasing. 
In the limiting case A —i- oo the deformation parameter e must vanish. Thus, 
the existence of an attractor solution 

qMr)^Br, BeR, (C.6) 

has been established. However, nothing can be said about the speed of the 
transition to the attractor solution and of course its soundness depends on the 
validity of the assumptions made to derive it. Finally, it has to be checked 
whether the constant B in (|C.6p has to scale to zero or not. 

Within this toy model it is even possible to provide a matter action explicitly: 

S{m) (q^Q + {dQ/dtf - U{Q)) drdt (C.7) 

liU{Q) is such that the sohxiion oi d'^Q / dt^ + U' {Q)+0{q^) = 0, Q{t), coincides 
with (1 — e) for all then the patching procedure above reproduces the exact 
dynamics of the whole system close to the asymptotics r = (which is the region 
relevant for the observer), because the equation of motion q + q{\ — Q) = 
from IjC.ip and (|C.7|I becomes equivalent to the one from ljC.3(l . However, it is 
emphasized that for the BH evaporation discussed in the rest of this paper no 
attempt is made to construct a corresponding matter action. 
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